Coulomb correlations do not fill the hole pockets in Nao.3Co02 
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There exists presently considerable debate over the question whether local Coulomb interactions 
can explain the absence of the small e'g Fermi surface hole pockets in photoemission studies of 
Nao.3Co02. By comparing dynamical mean field results for different single particle Hamiltonians 
and exact diagonalization as well as quantum Monte Carlo treatments, we show that, for realistic 
values of the Coulomb energy U and Hund exchange J, the e'g pockets can be slightly enhanced or 
reduced compared to band structure predictions, but they do not disappear. 



The Fermi surface of a material is one of its most fun- 
damental properties. Usually, it can be understood, at 
least qualitatively, within density functional theory. It 
came as a surprise, therefore, when several angle-resolved 
photoemission (ARPES)studies on the intercalated layer 
compound Nao.3Co02 [ll revealed a fundamentally dif- 
ferent shape of the Fermi surface than predicted by local 
density approximation (LDA) band theory 0. Recent 
bulk sensitive Shubnikov-de Haas data pi] also appear 
to be inconsistent with these calculations. On the other 
hand, the overall width of the Co 3c? bands in the ARPES 
data was found to be only moderately reduced compared 
to the LDA value. Essentially, the partially filled Co 
3d t2g bands should give rise to a large Og hole pocket 
centered around F, and six small hole pockets of e'g char- 
acter along the TK directions of the hexagonal Brillouin 
Zone. These e'g pockets have not yet been observed in 
experimental work. Their role for the superconducting 
hydrated phase of Nao.3Co02 is also a subject of intense 
investigations [3| . In view of the narrow width of the Co 
t2g bands {W « 1.5 eV), one possible source of the dis- 
crepancy between ARPES and band theory might be the 
effect of intra-3c? Co Coulomb interactions which, in prin- 
ciple, could enhance orbital polarization by leading to a 
charge transfer from ag to e'g subbands, and, eventually, 
to a shift of the e'g bands below the Fermi level. 

The influence of Coulomb interactions on the topol- 
ogy of the Fermi surface of Nao.3Co02 has been investi- 
gated by several groups, using various theoretical meth- 
ods and levels of approximation. Ishida et al. @ applied 
dynamical mean field theory (DMFT) [ej based on the 
multi-orbital quantum Monte Carlo (QMC) method, to- 
gether with a single-particle Hamiltonian derived from 
an accurate tight-binding fit of the t2g bands to the lin- 
earized augmented plane wave (LAPW) band structure. 
The result of this work was that, for Coulomb ener- 
gies U ~ 3.0 . . .3.5 eV and exchange J = C//4, the e'g 
hole pockets were slightly enlarged compared to the LDA 
Fermi surface, in contrast to the ARPES data. The width 
of the t2g bands, however, was found to be reduced to 
about 1 eV, in approximate agreement with ARPES. To 
avoid sign problems, only Ising-like exchange terms were 
included in the QMC calculation. It was also shown that 
an LDA-hU treatment can lead to enlarged or reduced 
e'g hole pockets, depending on whether U is smaller or 



larger than 5 J, respectively. 

At the same time, Zhou et al. [1] investigated this prob- 
lem within the Gutzwillcr approach in the large U , J — 
limit. Using a slightly different tight-binding fit to the 
LDA bands, these authors found that the e'g bands were 
shifted below Ep, and that the width of the t2g bands was 
strongly reduced from 1.5 eV to about 0.5 eV. Thus, while 
the Fermi surface appears to agree with the ARPES data, 
the band narrowing is much stronger than experimentally 
observed. Similar Gutzwillcr calculations in the U — *■ oo, 
J = limit were recently carried out by Shorikov et 
al. 0, with results similar to those of Ref. Q. Since 
the Gutzwillcr method replaces the frequency dependent 
complex self-energy by parameters providing orbital de- 
pending energy shifts and band narrowing, it represents 
an approximation to DMFT. Moreover, for [/ — > cx), com- 
plete orbital polarization is to be expected. Thus, for 
a meaningful comparison with ARPES data, it is im- 
portant to extend the Gutzwillcr approach to realistic 
Coulomb and exchange energies appropriate for Co. 

The influence of correlations on the electronic prop- 
erties of hydrated Nan35Co02 were also investigated by 
Landron and Lepetit [lOj within quantum chemical meth- 
ods for embedded CoOe and C02O10 clusters. The crys- 
tal field splitting between ag and e'g orbitals was found 
to be A = 315 meV, and the Coulomb and exchange en- 
ergies U = 4.1 eV and J = 0.28 cV. At present, it is 
not clear how these parameters, in particular, the large 
value of A, would be modified for larger clusters that are 
required to describe the electronic properties of the ex- 
tended system. Slave-boson mean field calculations by 
Bourgeois et al. [llj based on A = 315 meV and U ^ 00 
revealed a pure ag Fermi surface and a t2g band width of 
0.5 eV, similar to the results of Ref. [1]. 

To examine the role of Hund exchange contributions 
not included in the QMC/DMFT calculations, Per- 
roni et al. [l^l applied a new multi-band exact diag- 
onalization DMFT scheme to Nao.3Co02. This ap- 
proach does not suffer from sign problems and includes 
density-density contributions as well as spin-flip and pair- 
exchange terms. Also, larger values of U can be handled 
than via QMC. The result of this study was that in this 
material there is little difference between Hund and Ising 
exchange, and that, for U = 3... 5 eV and J — U/i, 
the e' pockets were slightly enlarged, in agreement with 
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the QMC/DMFT results. Also, the band narrowing was 
found to be roughly 30 %, consistent with the QMC treat- 
ment and with the ARPES data. 

Most recently, Marianetti et al. [l^ studied the prob- 
lem of the e'g hole pockets in Nao.3Co02 by applying a 
new continuous-time QMC/DMFT version that allows 
to reach larger U and lower temperatures. The single- 
particle Hamiltonian was the same as in Ref. [s^], ex- 
cept for the crystal field splitting A — Ea^ — E^'^ that 
shifts the ag bands up and the e'g bands down. For 
?7 = 3 . . . 5 eV and J = 0, reduced e'g pockets are found 
for A = — 10 meV, and fully suppressed pockets if A is 
increased to 50 . . . 100 meV. According to the authors, 
their results are "in agreement with Ref. ^] and in dis- 
agreement with Ref. 0". Since the DMFT calculations, 
however, were not done for the same input Hamiltonian 
and [/, J values as those in Refs. [1, the origin of the 
conflicting trends is presently unknown. 

The purpose of this work is to resolve this issue and to 
analyze the role of the single-particle Hamiltonian H{k) 
and Coulomb and exchange energies for the charge trans- 
fer between t2g bands. In particular, we show that the ED 
and QMC many-body calculations are in perfect agree- 
ment if identical input parameters are employed. On the 
other hand, the two different versions of H{k) used in 
Refs. [i[l3] and (below we refer to them as Hi and 

H2, respectively) give rise to a slight, but significant dif- 
ference in the variation of e'g occupancy with U: Whereas 
Hi yields decreasing e'g occupancy with increasing f/, 
H2 gives the opposite trend. We show that these dif- 
ferences are caused by the t2g crystal fields A contained 
Hi and i?2- The key point, however, is that, for realis- 
tic Coulomb and exchange energies, i.e., U ~ 3. ..5 eV 
and J w 0.72 eV the differences caused by H{k) are 
small and do not affect the controversy concerning the 
shape of the Fermi surface. Both versions of H{k) yield 
the result that, without an additional ag/ e'g crystal field 
splitting. Coulomb interactions do not eliminate the e'g 
hole pockets. The overall topology of the Fermi surface 
remains the same as predicted by LDA band theory. 

Fig. 1 shows the tight-binding fits to the partially occu- 
pied Co id t2g bands used in Refs. and [8, 13]. The 
total occupancy is 5.3. The ag and e'g occupancies (per 
spin band) are Ua ~ 0.80 and ng' « 0.925. Although 
both Hamiltonians give similar energy bands, they differ 
in a fundamental aspect: The predominant ag wave func- 
tion character of the lowest LAPW band along MK 
is correctly reproduced via Hi, resulting in a van Hove 
singularity in the ag density of states near —1.12 eV. In 
contrast, the two lowest H2 bands along MK cross, so 
that this singularity is shifted to —0.84 eV, implying a 
significant effective narrowing of this subband. The low- 
est H2 band at M has e'g character and gives only a weak 
step at —1.02 eV in the density of states. Since the in- 
fluence of Coulomb interactions is highly sensitive to the 
band width and the distribution of spectral weight within 
a band, these differences should affect also the correlation 
induced charge transfer between t2g bands. 
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FIG. 1: Upper panel: Tight-binding fits to LDA f2oband 
structure of Nao.3Co02 used in Refs. j^iO and Refs. [1, [Ill- 
Lower panels: comparison of ag and e'g density of states. Solid 
(red) curves: H\, dashed (blue) curves: H2, Ef = 0. 



The top of the H2 ag band at F is seen to exhibit a min- 
imum which is absent for Hi and leads to a pronounced 
peak in the density of states. This spectral weight is dis- 
tributed over a slightly wider energy range in the case 
of Hi . The minimum is caused by interlayer interactions 
which are absent in Hi as well as i?2- Moreover, the bulk 
Ug density of states does not exhibit a sharp peak in this 
region. Thus, the Hi ag density shown in Fig. 1 should 
be more appropriate than the one derived from H2. 

Despite these differences, near Ep both Hamiltonians 
yield nearly identical bands. The e'g bands extend less 
than 100 meV above Ep, and both models exhibit the 
ttg j e'g crossing along YK just below Ep. Note that these 
bands hybridize away from this symmetry direction, i.e., 
the crossing turns into an increasing hybridization gap as 
soon as the parallel momentum deviates from VK. 

We now discuss the correlation induced changes of the 
t2g bands of Nao.3Co02 as calculated within DMFT. We 
had previously demonstrated that, for the Hamiltonian 
Hi, the QMC and ED results of Refs. [1, [13] are in ex- 
cellent agreement and that both schemes yield reduced 
orbital polarization with increasing \J . Moreover, this 
trend was found to be insensitive to the choice of J. 
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We have applied the ED approach of Ref. to H2, in 
order to check its consistency with the QMC formahsm 
used in Ref. For [/ = 3 eV, J = 0, the subband 

self-energies as a function of Matsubara frequency were 
found to be in almost quantitative agreement. In view of 
the inevitable slight numerical differences between these 
fundamentally different DMFT approaches, the excellent 
agreement between the ED and QMC results is indeed 
remarkable. We also point out that both DMFT schemes 
take proper account of static and dynamical correlations. 

The unexpected result of this calculation is that, with 
H2 as input, both ED and QMC schemes yield enhanced 
orbital polarization: For U = 3 eV, J = 0, the subband 
occupancies are — 0.735, n^' ~ 0.957, compared to 
the LDA values Ua = 0.8, Ue' = 0.925. This charge 
transfer is opposite to the reduced orbital polarization 
obtained for Hi: — 0.825, Tie'^ = 0.91. 

Fig. 2 shows that similar systematic differences be- 
tween Hi and H2 are found at other values of U, J. To 
eliminate other sources of possible differences, all results 
are derived using the ED /DMFT approach of Ref. at 
T = 20 meV. Regardless of the choice of J, Hi leads to 
a reduction of function of U, whereas H2 yields 

increasing 71^1 . Thus, although the single-particle band 
structure ancf density of states derived from Hi and H2 
look qualitatively similar, these two Hamiltonians lead to 
a small, but significant difference in the variation of the 
subband occupancies with Coulomb energy. 

To analyze the origin of this unusual behavior we sim- 
plify the evaluation of the quasi-particle Green's function 

G{iLUn) = J2 [*^" + f^~ ^(^) ~ S(zu;„)]-i, (1) 

k 

where w„ — {2n+ 1)tt/ (3 are Matsubara frequencies, with 
(3 = l/ksT and temperature T. G, H and S are matrices 
in the t2g basis. Because of the planar hexagonal sym- 
metry, the diagonal elements of G are identical, and so 
are the off-diagonal elements. The same applies to E. In 
the Qg, e'g basis, these quantities become diagonal, with 
elements Ga = Gn + 2Gi2 and G^ = Gii — G12, and 
analogous expressions for S^, e- In this basis, the Green's 
functions can be approximately written as 

Gi{iuJn) = Jdujpi{uj)[iuJn + fl- LU - Tji{iuJn)]~'^ , (2) 

where Piiuj) are the Ug and e'g density of states compo- 
nents shown in Fig. 1. We have checked that Eq. (2) 
yields self-consistent DMFT solutions nearly identical to 
the ones derived from Eq. (1). Thus, the different solu- 
tions obtained for Hi and H2 are directly related to the 
different shapes of the respective density of states distri- 
butions. As is evident from Fig. 1, the e'g densities are 
quite similar for both Hamiltonians. Indeed, replacing 
one by the other does not alter the trends for the charge 
transfer shown in Fig. 2. It is clear, therefore, that the 
different shapes of the Og density of states are the source 
of the opposite orbital polarization found for Hi and H2 ■ 
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FIG. 2: Subband occupancies as a function of U for fixed 
J derived within ED/DMFT for T = 20 meV. Solid (red) 
curves: Hi, J = 0.7 eV; dotted (magenta) curves: Hi, J — 
0.9 eV; dashed (blue) curves: H2, J = 0.7 eV; long-dashed 
(green) curves: H2, J = 0.9 eV. For comparison, some results 
for J = are also shown (crosses). 



As pointed out above, the fact that the lowest H2 
bands cross along MK leads to an upward shift of the 
lowest Og van Hove singularity by about 0.3 eV. One can 
simulate this redistribution of spectral weight by reduc- 
ing the flg density of Hi in the range w < —0.7 eV and 
amplifying it in the region — 0.7eV< w < 0, such that the 
occupied weight remains 0.8. This deformation is suffi- 
cient to reverse the trend of ni{U) and give rise to a weak 
enhancement of orbital polarization. 

Evidently, the upward shift of spectral weight caused 
by the band crossing along MK implies a relative shift 
of 3d energy levels. In the case of Hi, the centroids of 
the Qg and e'g density of states are Ea^ = —0.624 eV and 
Ee' = -0.491 eV, where A = Ea ~ E^' = -133 mcV 
is the t2g crystal field splitting. In the tight-binding fit, 
A was varied along with the hopping parameters, in or- 
der to achieve the optimum representation of the LAPW 
bands throughout the Brillouin Zone @. Clearly, its 
value refiects the electronic structure of the extended 
system. In the case of H2, the splitting was chosen as 
A = — 10 meV, and only the hopping parameters were 
fitted @ . The Og and e'g centroids therefore nearly coin- 
cide: Ea = -0.489 eV and E^' = -0.479 eV. Note that, 
in both cases, Ea^ < Egi^ despite n^g < rie'^ . 

As a result of these different level splittings, correla- 
tions lead to an intriguing reversal of interorbital charge 
transfer: For Hi with A = —133 meV, the large Ug/e'g 
splitting is enhanced and gives rise to a gradual filling 
of the ttg band with increasing U. Since, at small U, 
the Qg occupancy is lower than the e'g occupancy, this 
charge transfer amounts to an initial reduction of orbital 
polarization. (At larger U, might become larger than 
Ue' , so that the same correlation induced e'g — > Og charge 
transfer eventually could turn into enhanced orbital po- 
larization.) In contrast, the small crystal field included 
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in iJ2, A = — 10 meV, is too weak to enforce a correla- 
tion induced downward shift of the Og band. Thus, the 
charge transfer is dominated by the larger occupancy, 
giving enhanced orbital polarization even at small U. 

Since Hi provides the more accurate fit to the LAPW 
bands, the correlation induced reduction of orbital po- 
larization in the range of reasonable values of U and J 
should be more realistic than the opposite trend obtained 
for H2. Thus, as argued in Refs. [1, correlations 
slightly enhance the e'g hole pockets of Nao.3Co02. We 
emphasize, however, that, according to Fig. 2, the oppo- 
site charge transfer obtained for H2 is not large enough 
to push the e'g bands below Ep. 

As pointed out in Ref. [13] , a large positive crystal field 
can enhance orbital polarization, so that, in combination 
with local Coulomb interactions, the e'g hole pockets dis- 
appear [1^. If we assume the ARPES data to be correct, 
the crucial question then concerns the physical origin of 
such a crystal field. Evidently, it is not related to on- 
site Coulomb interactions in the spirit of a single-site 
DMFT. Non-local effects stemming from the momentum 
dependence of the self-energy have not yet been explored 
and could be studied by using a cluster extension of the 
DMFT. Na disorder was shown to eliminate the pockets 
at large Na concentrations near x = 0.7 [IB, but is 
believed to be too weak to have a significant effect on 
the Fermi surface near x = 0.3. Surface effects which 
have played an important role in ARPES data on other 
perovskites, such as Ca2-2:Sra;Ru04 and Cai-ajSra-VOa, 
should also be investigated, in particular, the effect of 
Na induced states on the first layer. Moreover, possible 
structural distortions, such as intra-planar buckling, and 
their connection to the opening of the ag/e'g hybridiza- 



tion gap along TK should be explored. 

We finally mention that, as emphasized in Ref. [l^l, the 
e'g hole pockets are also important for the understanding 
of the heat capacity of Nao.3Co02. At present, the ex- 
perimental value, 7 12 ... 16 mJ/molCoK^, is difficult 
to reconcile with the LDA result, 14 mJ/molCoK^, 
and an effective mass enhancement of about 2, as derived 
within the DMFT studies discussed above. For a more 
detailed analysis of this quantity it might be necessary 
to relax the pinning condition implied by the single-site 
approximation and allow for non-local effects. 

In summary, we have resolved the puzzling discrepan- 
cies between DMFT results for the correlation induced 
dg/e'g charge transfer in Nao.3Co02. For identical in- 
put quantities, ED and QMC impurity treatments are 
in excellent agreement. Surprisingly, however, slight dif- 
ferences among the tight-binding Hamiltonians lead to 
increasing or decreasing orbital polarization. This rever- 
sal of subband occupancies as a function of U underlines 
the importance of using a high-quality single-particle ba- 
sis as input in the many-body calculation. In the present 
system, for realistic Coulomb and exchange energies, the 
differences resulting from these opposite trends are small, 
i.e.. Coulomb interactions do not fill the e'g hole pockets. 
The topology of the Fermi surface of Nao.3Co02 therefore 
remains the same as predicted by LDA band theory. 

We hope that these results encourage further experi- 
mental work to study in more detail the geometrical and 
electronic structure of this fascinating material. 

One of us (A.L.) likes to thank Chris Marianetti for 
extensive correspondence and for sending his QMC self- 
energy results, and Michelle Johannes and Igor Mazin for 
very useful discussions. 
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